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■ Abstract 



We compare gap equation predictions for the spontaneous breaking of global symmetries in 
supersymmetric Yang-Mills theory to nonperturbative results from holomorphic effective action 
techniques. In the theory without matter fields, both approaches describe the formation of a 
f — ■ gluino condensate. With Nf flavors of quark and squark fields, and with Nf below a certain 

critical value, the coupled gap equations have a solution for quark and gluino condensate for- 
' mation, corresponding to breaking of global symmetries and of supersymmetry. This appears 

to disagree with the newer nonperturbative techniques, but the reliability of gap equations in 
^ ■ this context and whether the solution represents the ground state remain unclear. 

p ■ 1 Introduction 



Spontaneous breaking of global symmetries in a gauge field theory sets in only when the gauge forces 

become strong. Accordingly, nonperturbative methods must be brought to bear, and the problem 

remains only partially understood. Approaches include lattice methods, semiclassical methods, 

and 1/N expansions. In addition, a simple approach based in Feynman graphs has been used with 

apparent success for many years. It comes in several closely related forms, including Schwinger- 

Dyson or gap equations, and "most attractive channel" analyses. These share the common idea 

that the forces responsible for symmetry breaking can be computed in a loop expansion. The lowest 

order contribution is typically one-particle exchange, with higher order contributions given by two 
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particle irreducible graphs jj], 0|. This can be applied to a most attractive channel force analysis or 
equivalently the kernel of a gap equation, or the CJT effective potential H for composite operators. 

Because the coupling must be relatively strong, there is no obvious small parameter in the 
expansion and this approach has been met with justifiable skepticism. Some evidence in favor of its 
utility, however, is provided by estimates of the second order contribution to the kernel 0, Bj. When 
the coupling is just strong enough to trigger the breaking, second order corrections are relatively 
small (less than 20%). With this check in mind, and in the absence of clearly superior alternatives, 
gap equations have been widely used in recent years in the study of dynamical electroweak symmetry 
breaking ||. Still, their validity remains open to question. 

Recent advances M in nonperturbative methods for supersymmetric gauge theories provide 
an independent framework in which to test the validity of the gap equation approach. The new 
methods lead to conclusions about spontaneous breaking of global symmetries in supersymmetric 
theories, which can be compared to gap equation results. This letter reports the results of such a 
comparison, first for the case of J\f = 1 supersymmetric Yang-Mills SU(N C ) theory with no matter 
fields (pure SQCD), and then for J\f = 1 SQCD including matter fields. 

2 Pure Supersymmetric QCD 

The Lagrangian for pure SQCD contains only gluon and gluino degrees of freedom (in Wess-Zumino 
gauge after eliminating auxiliary fields): 

Csym = -~ G a ^ v G^ va + J^G%& va + i X a ipX a + gauge fixing + ghosts, (1) 

where the canonically normalized Majorana spinor A describes the gluino field, 8 is the vacuum angle 
for the nonabelian gauge field G^ u , and iD^X = id^X — g [A^, X] . This theory has an anomalous 
U(1)r global symmetry, A — ► exp (ii?)A, broken by instanton effects to Z2N C - 

Various studies indicate the formation of a gluino condensate ]?], |8| , which spontaneously breaks^] 
the discrete Zin c symmetry to Z2. Nonperturbative holomorphic calculations derive this result 
starting from an exact Wilsonian effective superpotential for Nf = N c — 1, obtained |JSj, [K], 11] 



Note however a suggestion to the contrary 



from supersymmetric and global symmetry constraints, assuming the existence of a nonpertur- 
bative, supersymmetric regulator. Successively integrating out massive flavors results in an 
effective superpotential for Nf = 0. Differentiating with respect to the supersymmetric coupling 
t = iAk I g\ + 9 /2ir then provides an exact formula || for the gluino condensate, (A^A^) = 32-7r 2 A|. 
Here, the subscript h denotes a "holomorphic" normalization |l3| (1/g 2 , on the kinetic term) for 
a Weyl spinor, and = M exp [2irir / (3N C )] is defined in terms of a holomorphic gauge cou- 
pling g\ defined at an ultraviolet scale M in the dimensional reduction scheme For a canon- 



ical field normalization, rescaling A^ = g(M ) A requires the canonical coupling g to satisfy [14] 
8ir 2 /gl(M 2 ) = 8n 2 /g 2 (M 2 ) + N c lng 2 (M 2 ), due to the transformation's anomalous Jacobian p3| . 
We can accordingly express the holomorphic result as 

1 ,-„,„. 32vr 2 A o 32vr 2 A , 

5<**> = 7W) Ai = 7W> Al (2) 

(with canonically normalized, Majorana spinors). The one-loop (5 function solution for g diverges 
at the "confinement" scale Ai = Mexp [-87r 2 /(3iV c g 2 (M 2 ))]. 

Gap equation techniques, too, indicate a Z2Ar c -breaking condensate, signalled by a dynamical 
gluino mass. We write the gluino inverse propagator as A\(p 2 ) \jj> — Sa(p 2 )]* - > with the wave 
function factor A\(p 2 ) defined for renormalized fields so that A\(/i 2 ) = 1, and with the dynamical 
mass T,\(p 2 ) = A\(p 2 ) T,\(p 2 ). In a gauge with gluon propagator —i[g^ u + (£ — l)k^k u /k 2 ]/k 2 , the 
gap equation in Euclidean space takes the form 

~ , 2 , „ dl 2 , r dk 2 k 2 A X (M 2 ) Z\(k 2 ) 

where M and m are respectively the larger and smaller of k and p. (The angular integrations have 
been performed in four dimensions, after approximating — p] 2 ) by A\(Ai 2 ), and inserting an 

extra factor of A\(M 2 ) in anticipation of the linearization discussed below.) The structure of this 
equation is similar to that for the quark in ordinary QCD ||, ||], differing only in that the group 
factor = N c replaces C| und = (N 2 — 1)/{2N C ). We have neglected running of the coupling 
and of the gauge parameter ||, a crude approximation adequate to establish the existence of a 
symmetry-breaking solution in this theory and to set the stage for the presence of matter fields. 
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Here, the approximation requires an ultraviolet cutoff on the integral. When matter fields are 
included, an infrared fixed point will govern the transition and justify the neglect of running. 

A nontrivial solution to Eq. (||) requires the coupling to exceed a critical value, near which the 



dynamical mass vanishes continuously [15]. This allows the critical coupling to be determined by 
analyzing the equation in the regime p 3> T,\(p 2 ), where loop momenta k 3> T,\(k 2 ) dominate the 
integral and linearization in T,\ is a good approximation. We can then also neglect corrections to 
the massless renormalization group formula A\{p 2 ) ~ (// 2 /p 2 ) 7A , where the scaling exponent equals 
the gluino-field anomalous dimension 

7A = cfe I |^. (4) 

We find solutions by inserting a scaling form T,\(p 2 ) ~ Xa(/^ 2 ) (/J?/p 2 ) bx . The scaling ex- 
ponent b\, in leading approximation just (half) the gluino-mass anomalous dimension, obeys a 
^-independent equation to first order in g 2 : 

b x (i-b x ) = 3cf-^ + otf). (5) 

A solution requires the right hand side of this equation to be large enough for b\ to become complex 
|| , corresponding to the mass anomalous dimension reaching unity. The discrete global symmetry 
is predicted to break if the coupling exceeds the resulting critical minimum value, a cr = g 2 r /47r = 
ir/(3C^) = tt/(3N c ). Higher order corrections to the kernel of the gap equation may be small 
enough not to affect this conclusion qualitatively. In nonsuper symmetric QCD, the next order 
corrections are less than 20% [||, [|]. 

In reality, the coupling runs; this effect can be incorporated in a WKB-like approximation || 
which predicts nonzero solutions at scales where the coupling exceeds the above critical value. Since 
the theory is asymptotically free and presumably confining, that will occur at momenta of order 
the confinement scale.0 The symmetry is thus predicted to break, and the confinement scale sets 

the size of T,\ via the full nonlinear gap equation. Inserting this result into the operator product 

2 The Wilsonian effective coupling has been argued to be controlled by an exact j3 function Jl^] for Nf = 0: 
dg/dln/i = — (37V c /167r 2 ) <? 3 /(l — N c g 2 /(8ir 2 )). Well before reaching the singularity in this (3 function, the running 
coupling will reach the above g CT . The f3 function's expansion parameter there is N c a CI / (2tt) = 1/6. 
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formula for the gluino condensate leads to the estimate 

The last rough estimate, neglecting logarithmic factors, arises from dimensional analysis with a 
nonzero solution to the gap equation, identifying the confinement scale with Ai in Eq. (0). For 
large N c , we may write this estimate in the form 

which has the same N c scaling (with fixed N c g 2 (M 2 )) as Eq. (|2|). 

We note that a complete gap-equation prediction of symmetry breaking would require comparing 
the broken solution to the = solution. Establishing the former to be the ground state would 
require settling currently unresolved issues of gauge dependence, resulting from truncating the 
effective action or the infinite set of coupled Schwinger-Dyson equations. (Even in QCD, gap 
equation techniques face the same question.) Assuming the broken solution to be the ground state 
(with zero energy for unbroken supersymmetry) , the gap equation and holomorphic approaches 
both predict that the discrete Z2N C symmetry spontaneously breaks to Z2 via a gluino condensate. 

3 Supersymmetric QCD with Nf Flavors 

Including quark supermultiplets with Nf flavors provides more scope for comparing the two ap- 
proaches. Again using component notation with auxiliary fields eliminated, we add Nf quark 
flavors to the Lagrangian, as Dirac spinors in the SU(N C ) fundamental representation; their left- 
and right-handed components are respectively associated with scalar superpartners, <f> and 

Csqcd = C SYM + H'ip-mo)^ + \D^\ 2 + \D^\ 2 - mg(0*0 + #*) 

- y (V T a (/> - ^T a 0*) 2 - igV2 (j)*T a \ a P L i) + 4>T a \ a Pf{ip - h.c.) , (8) 

where iD^ = id^— gA^T 0, . In the limit of vanishing bare masses mo — > 0, the global symmetries are: 
a discrete parity symmetry, (ipL,<f>) <-» (ipR,<f>*)\ chiral SU(Nf) x SU(Nf) with quarks and squarks 
rotated together; baryon number U(1)b for quarks and squarks; and the anomaly-free subgroup 
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U(l)x of the gluino/squark U(1)r and the quark/squark axial U(1)a- Although these symmetries 
forbid perturbative mass generation, nonperturbative condensates could break the symmetries, 
generating fermion and scalar dynamical masses and mixing the left- and right- scalars. The scalar 
mixing would also require nonvanishing dynamical gluino mass, since U(l)x combines gluino and 
axial rotations. 

An important step in the nonperturbative study of supersymmetric theories was Witten's ob- 



servation [17] that in a massive theory the difference between the number of bosonic and fermionic 
zero-energy states is conserved^] (the Witten index), with a nonzero value (N c for massive SQCD) 
implying unbroken super symmetry. In that case, supersymmetric Ward identities for Nf > leave 
fermion condensates proportional to the bare mass. Specifically, for the gluino condensate the 
Konishi anomaly [19] gives \ (AA) = (g 2 / 16tt 2 ) mo (<j)<j))- For the quark condensate [pCfl , each flavor 



i satisfies {ipiipi} = mo {4>*<j)i + &</>*)■ The massless theory thus forbids fermion condensates, as- 
suming a well-defined, nonvanishing index at mo = 0, and assuming squark bilinear expectations 
are not too singular when mo — ► 0. 

The holomorphic effective action approach || provides a more complete description of SQCD 
with Nj flavors. Assuming the existence of a nonperturbative, supersymmetric regulator, it leads 
to a self-consistent picture for various values of Nf. Many physically distinct but (if Nf > N c ) 
degenerate vacua form a moduli space, on which generically the chiral symmetry breaks, but at 
special points (when expectation values vanish) is preserved. 

• For Nf > 3N C , there is no asymptotic freedom, and the quarks, squarks, gluons and gluinos 
become noninteracting at large distances. 

• For |iV c < Nf < 3N C , the theory is asymptotically free and the coupling runs to a "supercon- 
formal" infrared fixed point. There is no confinement in this "nonabelian Coulomb phase", and 
the global symmetries remain unbroken at the origin of moduli space (where the squark vacuum 
expectation values vanish). The spectrum is described by interacting, massless quarks, squarks, 
gluons, and gluinos. 

3 For the exactly massless theory, zero modes can come in from infinity and leave the index ill-defined. This is 
avoidable Jx8| if Nf is an integer multiple of N c . 
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• For N c + 2 < Nf < |iV C! there is no confinement, but the gauge theory is strongly interacting in 
the infrared. The spectrum is best described in terms of massless, composite mesons and baryons of 
a local, infrared-free, dual gauge theory. At the origin of moduli space the theory leaves unbroken 
all the original global symmetries. 

• For Nf = N c or Nf = N c + 1, the theory confines and the spectrum consists of massless compos- 
ite particles corresponding to meson and baryon fields composed of the original matter fields. For 
Nf = N c , either U(1) B or chiral SU(N f ) x SU(Nf) must break. For Nf = N c + 1, there can be 
confinement without chiral symmetry breaking. 

• For 1 < Nf < N c , a nonzero superpotential |11] lifts the degeneracy and fixes the ground state 



at arbitrarily large scalar expectation values (for mo — * 0). 

An important feature of these results is that for Nf > N c , the global symmetries associated with 
massless quarks, squarks and gluinos can all remain unbroken. 

We next use the gap equation framework to study the SQCD theory with matter. A previous 



investigation \21\ concentrated on the quark and squark dynamical masses, neglecting the gluino 
condensate by assuming the matter fields to occupy a high dimensional representation. Here we 
will retain all condensates, to obtain a set of coupled gap equations for the quarks, gluinos, and 
squarks. The most attractive interactions occur in channels that preserve local SU(N C ) and parity, 
so we will not consider color-breaking condensates, or differing dynamical masses for the scalars <f> 
and 4> (Eq. |||). We will also exclude nonzero vacuum values for the squark fields. Supersymmetry 
requires equal bare masses for the quarks and squarks, which we take to vanish; but the dynamical 
masses can differ, since the component notation does not maintain manifest supersymmetry in 
off-shell Green functions. 

For small enough coupling, the running will be determined by the (scheme-independent) two- 



loop (3 function |22] 



(3N C -N f ) - *£—(3Nt-2N c N f + N f /N c ) + ••• , (9) 



din// 16vr 2 v J ' (16tt 2 ) 2 

which displays an infrared fixed point when 3A r c /(2 — N^ 2 ) < Nf < 3N C . It occurs at 

3N C -N f 

a * ~ 2?r - 3AT C 2 _ Nf/Nc > ( 10 ) 



a small value for Nf sufficiently close to 3N C . The coupling approaches this fixed point at scales 
below some intrinsic scale A governing the solution to Eq. (||). If the fixed point coupling exceeds 
a critical value, then the gap equations will exhibit nontrivial solutions. Since the dynamical mass 
vanishes continuously there,^ the transition will be governed by the fixed point, and the running 
of the coupling may be neglected to first approximation. Furthermore, the gap equations may be 
linearized in the neighborhood of the transition. 




<\At@w = vKtvfv + W — O — 




H 



Figure 1: Gap equations for the gluino and for the quark masses. Solid lines are quarks, dashed 
lines are squarks, wiggled lines are gluons, and wiggled-on-solid lines are gluinos. Blobs symbolize 
the full propagators and vertex structures. The scalar exchange graphs give integrands proportional 
to higher powers of dynamical mass, and are dropped in the linearized approximation. 

The gluino and quark gap equations are shown in Fig. 1. To lowest order in the kernel and 
in the linearized approximation relevant to the study of the transition and the determination of 
the critical coupling, they each receive a contribution only from the first graph (gluon emission 
and reabsorption) . The gluino gap equation is simply the linearized version of Eq. (|3|), and the 
corresponding quark equation takes the form ||| 

^( P )-gC 2 (3 + J i^WA^j^^' (U) 

where g 2 is now the value of the coupling at the infrared fixed point. For both, the intrinsic scale 
A is effectively an ultraviolet cutoff in this approximation: the functions S(p 2 ) steepen there from 
near-critical behavior ~ 1/p, to asymptotic behavior ~ 1/p 2 - The gauge parameter is also in general 

a function of momentum, governed by a renormalization group equation. When g 2 approaches its 

4 Ref. |l5| argued that although the dynamical mass vanishes continuously, there are no light degrees of freedom 
in the symmetric phase and therefore the transition is not, strictly speaking, of second order. 
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infrared fixed point, the gauge parameter has one as well. We denote this fixed-point value simply 
by £, and in this approximation take it outside the integral. The wave function factors A(p 2 ) are 
computed in the massless theory, dropping corrections subleading in mass. They are defined to be 
unity at the renormalization scale /x 2 . 

For the scalars, we obtain a gap equation for the diagonal mass E 2 (p 2 ), appearing in both 
(^Wf(-p)) and $*(p)0(-p)>: 



P 



\A^)-Z^) + ^(p 2 )S 2 (p 2 > 



(12) 



a 



fund 



dk 2 

16vr 2 



M k2 ) 



k 2 



M 2 > k 2 + E|(A; 2 ) 



+ (3 + MM 



^MM 2 ) 

AJk 2 ) 



MM 2 )MM- 



A^{k 2 ) k 2 + T 2 {k 2 ) 



{k 2 + ^r(k 2 - P 2 ) - Z 2 (M 2 ) • [1 + sign (k 2 - p 2 )]} 



where the wave function renormalization constant implements A ( f ) (fi 2 ) = 1, and A g (k 2 ) repre- 
sents the gluon wave function factor .0 The dynamical mixing mass Y? x {p 2 ) in {4 > {p) ( t ) {~P)) satisfies 



^fund 
°2 



dk 2 g 2 



16tt 2 M< 



A<f,(m 2 



Z 2 x (k 2 ){(Z + l)M 2 + (3-0m 2 } 



(k 2 + £ 2 (k 2 ) 



+4 M^m ^ ^ M k 2 + s|(fe 2 ) 



(13) 



In Eqs. (12) and (13) we have dropped some contributions from diagrams with higher powers of 
mass insertions, as well as higher order contributions from diagonalizing the scalar propagator. 
Supersymmetry ensures cancellation of the quadratic divergences in the diagonal mass equation 
fll2]). After that cancellation, both equations can be fully linearized for momenta p large compared 
to E. 



Then up to mass-suppressed corrections, A r p(p ) is fixed by Eqn. (12), while A\(p ) « (// jp 



and A 1 p(p 2 ) (fJ. 2 /p 2 )^ , where 



7A 

7v 



fund J 



16vr 2 



(£ + l)C£ 



fund 



(14) 



5 Here, g and £ are in general functions of M 2 . As in the fermion equations, this dependence can be neglected if 
the integrals are dominated by an infrared fixed point. 
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Here Cf un d = 1/2, and g 2 and £ are the fixed point values described above. Since the critical 
couplings depend on the wave function factors, the gluino critical coupling will differ from its value 
in the matter-free theory, even though in both cases gluon emission and reabsorption is the only 
leading contribution to the gap equation. 

We focus initially on the quark and gluino equations. As described for the matterless theory, 
substituting the wave function factors of Eq. into the gap equations, and setting the mass 
anomalous dimensions to unity, leads to critical couplings independent of £: 

ZCf-2N f C hmA 3N C -N f > 
~ 2C^ d ~ N c -1/N c ■ [ 1 

The choice of Nf determines whether the infrared coupling a* in Eq. (|To|) achieves these critical 
values. As Nf is reduced from 3iV c , a* increases from zero and first exceeds the quark critical 
coupling, a cr ^, when Nf/N c = (f ) il/q/^l v- a ~ 2.25. For this value of Nf/N c , a* does not 
yet reach the gluino critical coupling, leaving quark condensation to play the primary role here. 
Assuming validity of the approximations above, we can now outline the phase structure of SQCD 
as described by gap equations: 

• For Nf > 3N C there is no asymptotic freedom, in agreement with the holomorphic description. 

• For 2.25N C < Nf < 3N C there is an infrared fixed point at which the two-loop (3 function vanishes. 
For Nf close enough to 3-/V c , the fixed point coupling is small, keeping the coupling well below the 
critical values of Eq. (|l5"l). For any Nf above 2.25-/V c , the infrared coupling remains below both 
critical couplings, leaving the quark and gluino symmetries unbroken. Confinement does not set 
in and the theory remains in the nonabelian Coulomb phase, in agreement with the holomorphic 
prediction. 

• For Nf < 2.25A r c , the infrared coupling (in the two-loop approximation) exceeds the quark critical 
coupling. If higher order effects may be neglected, the gap equation then indicates that quarks 
develop a dynamical mass which vanishes continuously at the transition. The quarks decouple from 
the j3 function at this mass scale, eliminating the fixed point and causing g to run to the gluinos' 
(larger) critical coupling. Gluino and scalar dynamical masses are then generated at essentially the 
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same scale as the quark mass. Below these mass scales, decoupling allows g to continue running to 
confining values. The global chiral and U(l)x symmetries are thus predicted to break in this Nf 
range. The supersymmetric Ward identity then requires that supersymmetry breaks as well. 
Holomorphic methods, by contrast, predict the nonabelian Coulomb phase to remain down to 
Nf = |iV c , and allow preservation of the global symmetries without confinement even down to 
N f = N C + 1. 

These patterns also control the scalar gap equation solutions, which depend upon the fermion 
solutions. The inhomogeneous scalar equations do not admit vanishing solutions when the quark 



and gluino dynamical masses are nonzero [21]. In that case, the scalar solutions should be of 
comparable scales. The connection between scalar and fermion mass solutions can be understood 
from the global symmetries: a nonzero off-diagonal mass T, x (p 2 ) breaks both U(l)x and SU(Nf) x 
SU(Nf), and thus requires corresponding gluino and quark masses. The diagonal mass S^(p 2 ) 
breaks no chiral or U(l) symmetries, but is not forced to vanish unless the quark mass does and 
supersymmetry is unbroken. 

Checking that these symmetry breaking solutions of the gap equation describe the ground state 
would require showing that the resulting vacuum energy is lower than that of the chirally symmetric 
solution. (For this to be the case the symmetric solution, as well as the broken one, would need to 



break supersymmetry.) For SQCD with matter, suggestive calculations [21, 23] do indicate lower 
energy for the symmetry breaking vacuum, although these are gauge dependent and neglect squark 
and/or gluino dynamical masses. 

4 Conclusion 



In pure SQCD without matter fields, various approaches indicate that the discrete global symmetry 
%2N C breaks spontaneously to Z<i through the formation of a gluino condensate. The theory confines 
and supersymmetry presumably remains unbroken. In the gap equation, the lowest order kernel, 
single gluon exchange, provides an attractive force capable of condensate formation. Whether the 
broken solution represents the ground state of the theory, however, remains unclear. Assuming 
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this to be the case, the gap equation and holomorphic effective action techniques lead to the same 
conclusion. 

For SQCD with matter multiplets, gluon exchange alone again gives the leading contribution to 
the linearized gap equation governing the transition. At this order, the gap equation (with gauge 
coupling governed by the two-loop (5 function) indicates that spontaneous breaking of the global 
symmetries occurs for Nf below a critical value near 2.25N C . This transitional value corresponds to 
the infrared fixed point coupling exceeding the critical strength necessary for quark condensation. 
If we instead determine the fixed point coupling from the three-loop (3 function (in the DR scheme) 
[22], the corresponding critical Nf ~ 2.08N C , only an 8% shift in this renormalization scheme. 



Examining the size of higher order terms in the kernel of the quark gap equation, corresponding 
to the mass anomalous dimension at two loops, would give another check on the validity of these 
approximations. In nonsupersymmetric theories, such kernel corrections are of order 20% Q. 
With these caveats, for Nf just below the critical value the quark mass is nonzero but small, 
justifying the linearization. 

At scales below their predicted mass the quarks decouple, eliminating the infrared fixed point 
so that the coupling increases to the gluino critical value at essentially the same scale. The gluinos, 
too, thus condense at the quark mass scale. Squark masses are induced by the fermion masses, and 
are also of the same scale. The supersymmetric Ward identity then predicts that supersymmetry 
also breaks in this regime. An effective theory with Goldstone bosons and Goldstinos emerges at 
low energies. All this is in contrast with the picture emerging from holomorphic effective action 
techniques, where the nonabelian Coulomb phase is argued to persist down to Nf = ^N c and the 
global symmetries to remain unbroken (at the origin of moduli space) down to Nf = N c + 1. 

Of course, as already stressed, gap equation solutions in general have not been shown to corre- 
spond to ground states. (Supersymmetric constraints on ground state energies may be important 
in this connection.) Settling the question involves still-unresolved issues of gauge invariance and 
interpretation of the effective potential. Furthermore, truncating kernels at one loop, and the f3 
function at two loops, seems quantitatively less reliable than in nonsupersymmetric theories. The 
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kernel truncation, together with the use of component notation (Wess-Zumino gauge), may also in 
effect explicitly break supersymmetry. The latter possibility can be checked by analyzing a set of 
gap equations in the manifestly super symmetric superspace formalism p3]. A recent study of this 



problem [25] indicates that symmetry-breaking solutions could exist, although infrared divergences 



in the formalism have so far obscured the analysis. 

The possible discrepancy between gap equation analyses and results based on holomorphic 
effective action techniques and index theorems is noteworthy and invites further study. Efforts to 
reconcile them may help deepen our insight into the behavior of strongly coupled gauge theories. 
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